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Abstract. We describe a system for long-distance distribution of quantum en- 
tanglement, in which coherent light with large average photon number interacts 
dispersively with single, far-detuned atoms or semiconductor impurities in optical 
cavities. Entanglement is heralded by homodyne detection using a second bright 
light pulse for phase reference. The use of bright pulses leads to a high success 
probability for the generation of entanglement, at the cost of a lower initial fidelity. 
This fidelity may be boosted by entanglement purification techniques, implemented 
with the same physical resources. The need for more purification steps is well 
compensated for by the increased probability of success when compared to heralded 
entanglement schemes using single photons or weak coherent pulses with realistic 
detectors. The principle cause of the lower initial fidelity is fiber loss; however, 
spontaneous decay and cavity losses during the dispersive atom/cavity interactions can 
also impair performance. We show that these effects may be minimized for emitter- 
cavity systems in the weak-coupling regime as long as the resonant Purcell factor is 
larger than one, the cavity is over-coupled, and the optical pulses are sufficiently long. 
We support this claim with numerical, semiclassical calculations using parameters for 
three realistic systems: optically bright donor-bound impurities such as 19 F:ZnSe with 
a moderate-Q microcavity, the optically dim 31 P:Si system with a high-Q microcavity, 
and trapped ions in large but very high-Q cavities. This is a preprint. Please 
consult published version of paper freely available as New. J. Phys. 8, 184 (2006), 
at http -.//stacks .iop. org/1367-2630/8/184 ■ 
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Figure 1. The three-state level structure of the atom or impurity complex. The 
two long-lived ground states |0) and |1) form the qubit and feature a large energy 
separation HA. Only the transition between state |1) and state |e) is optically active. 
The optical pulse is detuned from this transition by uj p and the cavity mode is detuned 
by wo- 



1. Introduction 

The interaction of an intense, off-resonant optical pulse with a single atom in a cavity 
has been the subject of a number of experiments pQ. Recently, such interactions have 
become important for non-destructive measurement of atoms in weak traps. In such 
systems, spontaneous emission would evict the atom from the trap, so maximizing the 
detectable phase shift from an off-resonant interaction while minimizing absorption has 
been the subject of several recent studies pOElll]- 

If phase shifts large enough to detect the presence of a single atom with negligible 
absorption are indeed available, then these off-resonant interactions can form the basis 
of a robust means of distributing entanglement between atoms or impurities in distant 
cavities [Sj. To understand how, consider a qubit formed from the two ground states 
of a A-type system, as in figure ^ in which only one of the ground states interacts 
with the light. In this system, an absorption- free dispersive interaction is sufficient to 
measure the state of this qubit. However, if the light pulse is sent to a second cavity 
with a second qubit and only measured afterwards, and if the measurement shows a 
phase shift corresponding to one and only one qubit in the optically-active state, the 
lack of information as to which qubit was in which state may post-select an entangled 
Bell state. Such entanglement forms the basis of a quantum repeater; it is a "hybrid" 
system because it relies on both the discrete states of the individual A-type emitter and 
the continuous quantum variable of the coherent light amplitude. 

Realistically, the fidelity of this operation is reduced by any leaked "which-path" 
information. This loss may be provided by light leaking from the fiber between the 
cavities, and it may be provided by a small probability for atomic absorption during the 
supposedly dispersive interaction. The former effect is unavoidable, and ultimately 
limits the fidelity of the post-selected entangled state to a value much lower than 
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the theoretical expectations for most other proposals for CQED-based entanglement 
distribution [HI El El 13 113 IH] ■ However, in this case the light used is bright and readily 
available from a laser source, the detection can be done with high efficiency, and the 
critical phase information can be stabilized by sending a trailing reference pulse down 
the same optical fiber. For these reasons, the expected fidelity, though low, is realistic 
with existing fiber-based technology, and the probability of successful post-selection of 
the desired state is very large, leading to a very fast initial entanglement distribution. If 
we then add a suitable protocol for nested entanglement purification and entanglement 
swapping the low fidelity may be fairly quickly improved to near-unity values. 

The generation of entanglement in our scheme is probabilistic but heralded, like 
a number of previous proposals. Deterministic CQED-based entanglement distribution 
schemes exist [TTj , but such schemes put more challenging constraints on optical 
cavities and optical pulses. Duan et al |H] proposed a probabilistic scheme using Raman 
scattering from atomic ensembles. This scheme has recently seen an experimental 
implementation jTH]- Integration of this system into a larger repeater architecture may 
be difficult; proposals operating on the principle of state-selective scattering from single 
emitters in cavities may be better suited for realistic implementations. Childress et al [0] 
have described a scheme in which such emitters reside in each path of an interferometer. 
Duan and Kimble have proposed the use of the 7r-phase shift incurred upon reflection 
of a single photon from a loaded cavity as the basis for quantum logic ^3] , and a similar 
principle has been proposed for a quantum repeater design by Waks and Vuckovic [TU] . 
with an architecture consistent with photonic-crystal microcavities. In common to all of 
these schemes is a reliance on single photon (or sub-photon coherent state) transmission 
between the distant qubits. Although the heralded entanglement may have high initial 
fidelity, the high probability for channel loss means that the frequency of successful 
entangling events decreases exponentially with the distance between repeater stations. 
Since the rate of attempts is inevitably limited by the classical communication time 
between the distant stations, the use of single photons greatly reduces the speed of 
entanglement generation. The present scheme's crucial difference is the use of bright 
pulses which assure successful post-selection for about 36% of the pulses sent down the 
channel. The trade-off for this increased speed is the limited fidelity of the post-selected 
entanglement, requiring additional entanglement purification. 

For entanglement purification and swapping, local quantum logic is necessary. This 
step is especially challenging in proposals well suited to long-distance entanglement 
distribution, such as systems based on atomic ensembles 0. If large and nearly 
absorption-free phase shifts are available, however, the same off-resonant interactions 
used for entanglement distribution may be used for a controlled-sign gate based on the 
geometrical phase j!5j . This gate is measurement-free and completely deterministic, 
allowing rapid entanglement purification and swapping. However, it is far less robust to 
optical loss than the entanglement distribution scheme. 

The final entanglement fidelity will depend on the fiber loss, certainly, but also on 
the amount of absorption in the dispersive interaction between the optical pulse and the 
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emitter /cavity system. In this work we will show the physical criteria required for the 
desired qubit coherence to be well preserved even in the presence of this small but finite 
absorption. We will focus on systems in the intermediate coupling regime, or "bad cavity 
limit," where the light-matter coupling is smaller than the rate of light-leakage from the 
cavity, but strong enough to substantially modify the rate of spontaneous emission. This 
regime is relevant for practical, homogeneous emitters in solid-state cavity systems. 

We will begin in section [21 by considering how initial probability of success and 
fidelity affect the final rate of quantum communication; this discussion will show the 
potential advantage of the current proposal. Then section El will analyze the procedure 
for long-distance entanglement distribution and local quantum logic, in more detail than 
the sketch above or our previous treatment [5 . For this treatment the CQED interaction 
will be idealized. Once we have established the important figures of merit, we will 
estimate the ability for real CQED systems to implement this proposal. Our methods of 
analysis are explained in section HJ with results presented for three different regimes of 
operation in sectional 

2. Motivation: Final Rate of Communication after Entanglement 
Purification and Swapping 

The means of distributing long-distance entanglement sketched in the introduction 
differs from most other proposals for heralded entanglement primarily in the fact that 
the probability of successful initial entanglement generation is high (~ 36%), while the 
initial fidelity is low (~ 77%). Before analyzing the origin of these numbers, let us 
address the question of whether such a scheme is useful. Obviously, an increase in the 
probability of success will increase the final rate of communication. However, the need 
for more entanglement purification will also slow down the final rate of communication. 
The question we address in this section is how the probability of success and fidelity affect 
the rate of final long-distance entanglement in a complete quantum repeater architecture. 

2.1. Entanglement Purification and Swapping Protocol 

The answer to this question depends on the protocols used for entanglement purification 
and entanglement swapping. Three such protocols for entanglement purification were 
presented in This important work analyzed the efficiency of "nested purification" 
schemes, in which imperfect distance-doubling entanglement swapping procedures are 
followed by entanglement purification. These schemes consider N + 1 repeater stations 
(including end- stations), where N = £/£, for total distance C subdivided into distances 
I between adjacent stations. Very fast schemes ("A" and "B") were considered in which 
hundreds of logically connected qubits are present in each station, as well as a scheme 
( "C" ) in which as few as two qubits and at most 2 log 2 N are present in the stations. 
More recent work j5] has shown that it can be sufficient to have only two qubits in every 
station. In these protocols the initial fidelity of entanglement distribution is considered 
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to be quite high, so that purification is only used to correct for fidelity degradation and 
gate error during entanglement swapping. 

Schemes using a minimal number of qubits are inevitably slow, since the initial 
entanglement generation and purification protocols are probabilistic, and with a minimal 
number of qubits, the entire protocol must begin from the start after each failure. A slow 
protocol cannot be remedied by arbitrarily speeding up the entanglement generation and 
purification procedures, since these are inevitably limited by the time it takes to transfer 
information between adjacent stations. The rate of final entanglement generation does 
speed up considerably if an ensemble of qubits is present in each station, as in schemes 
A and B of ^2]- However, these schemes present an extreme case where arbitrarily 
many qubits are allowed in order to exponentially speed the purification process (the 
number needed grows polynomially with distance). As logically coupled qubits present 
an expensive resource, such schemes may be too expensive. 

We consider instead a scheme similar to scheme B of J2| > still using an ensemble of 
qubits, but we restrict the size of that ensemble as much as possible. From scheme 
C of [12], we know that at least 21og 2 iV qubits are needed to allow the distance- 
doubling nested purification scheme to proceed in parallel to entanglement generation. 
If the initial fidelity begins low, so that entanglement of qubits in adjacent stations 
requires initial purification, two more qubits are needed in each section. Therefore at 
least 2 + 2 log 2 iV qubits are needed. To speed the protocol without increasing qubit 
overhead more than this, we consider putting 2 + 2 log 2 iV qubits in each of the repeater 
stations (including the endpoint receiver and sender). Initial entanglement is generated 
in parallel between the 1 + log 2 N "send" qubits and the 1 + log 2 N "receive" qubits 
in each station. This parallel operation significantly improves the speed of the initial 
entanglement generation, and allows simultaneous generation of new entangled pairs 
while long-distance pairs wait for purification. 

With this number of qubits assumed, the protocol followed is straightforward. Each 
station purifies a certain number of steps as successful entanglement post-selections 
occur. Once the prescribed number of purification steps are done for both a send 
and a receive qubit in the same repeater station, Bell-state analysis is performed for 
entanglement swapping, assuring that each such process doubles the distance over 
which the qubits are entangled. After purification or swapping operations, each station 
immediately begins sending and receiving new optical pulses to and from any available 
qubit for further entanglement generation. The specific purification procedure we 
consider is similar to that used in scheme B of the recurrence protocol originally 
presented in ^B]- One slight modification which will be important for the specific 
entanglement scheme we discuss in later sections is that in each purification step, the 
local rotations performed before the controlled-NOT gates are chosen case-by-case to 
optimize the subsequent fidelity. This operation requires no additional complexity in 
the implementation as long as the form of the noise present in the system is known in 
advance. 
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2.2. Probability of Success vs. Initial Fidelity 

With such a protocol established, we may now begin to compare the interplay between 
initial probability of success and initial fidelity. We numerically simulate the protocol 
described above for N + 1 = 129 stations, assuming that all processes are limited 
only by a communication time of 50 /is between adjacent stations, corresponding to 
about i = 10 km station-to-station distance and a total distance of 1280 km. Once 
an entangled pair is generated at 1280 km, a quantum bit is teleported across that 
distance. When the teleported bit arrives at the last repeater station, the simulation 
marks the time, resets the qubit participating in the teleportation, and continues to 
work on generating more entanglement. It takes some time for the first long-distance 
pair to be generated; it takes less time for subsequent pairs to be generated as there is 
already some entanglement in the system. The differences in arrival times of between 
5 and 30 bits are averaged, and from this we obtain an average rate of long-distance 
quantum communication. 

This rate is very fast in the unrealistic case that the initial fidelity is high and no 
gate errors are present. A more realistic rate must consider the possibility of imperfect 
entanglement and faulty gates. In this section we model error in a general way: we use 
a simple white-noise model for a two-qubit density operator, defined by 



where 1 is the two-qubit identity matrix. 

Figure |2 shows results of such simulations, in all cases with a number of purification 
steps at each distance chosen to yield a final communication fidelity of about 95%. We 
have considered two cases in generating this plot. The first models proposals based 
on single photon detection [HI IHl HD] , where the initial fidelity is quite high. We have 
chosen that the initial state is a desired Bell state with white noise added at a level 
of e = 5%, resulting in an initial fidelity of 96%. This number is comparable to 
theoretical expectations, although it is much lower than fidelities observed in existing 
experiments jTSj. In schemes such as these, the probability of success is limited by 
channel loss, the total efficiency of the photon detection, the scattering probabilities 
from the qubits to be entangled, and a factor of 1/2 from the post-selection probability. 
Taken together, these factors can easily result in a realistic probability of success less 
than 0.001. In contrast, we have also considered a model in which the initial fidelity 
is quite low, using e = 30% resulting in 77% fidelity, a number we will motivate in the 
next section. At equal probability-of-success, such a scheme is of course much slower. 
However, we will see that it allows a realistic probability of success on the order of 36%, 
which is far out-of-reach of proposals based on single photons. Figure |2] shows us that if 
we compare to probabilistic single-photon based schemes with probability of successful 
entanglement generation less than about 0.5%, the much higher probability of success 
of 36% present in our proposal results in a higher final communication rate, despite the 
need for additional purification. 

The white-noise model used for entanglement generation and gate error in this 
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Figure 2. Final communication rate for a full nested purification protocol assuming 
an abstract white-noise process for the initial state generation and local gate error, 
with numbers of purification steps chosen to target a final fidelity of about 95%. The 
squares and circles, corresponding to an initial fidelity of 96%, show a final rate which 
varies slightly sublinearly with success probability (oc p - 93 ). The triangles correspond 
to an initial fidelity of 77% and show a slightly superlinear dependence (oc P 12 ). 



analysis leads to slower purification rates than what may be possible in actual practice. 
This model is only used here to isolate the issue of probability-of-success and fidelity 
from other issues. In section IH.5[ we will revisit this swapping and purification protocol 
using more realistic error models for our proposal, and we will see that the expected 
communication rates are somewhat faster than the results in figure El 



3. Quantum Information Processing with Dispersive Light-Matter 
Interactions 

In this section, we analyze the basic means of generating entanglement using dispersive 
CQED interactions with bright coherent light, including the dominant error modes. 



3.1. Idealized System 

The basic qubit in this scheme is formed by the two lower states of a three-state A- 
system, as shown in figure HJ The two metastable qubit states are labelled |0) and |1). 
We define the Pauli-Z operator for the qubit as 

Z=|0)(0|-|1)(1|. (2) 

Coherent transitions (rotations) between these two states are presumed to be possible 
through methods we will not discuss here, such as stimulated adiabatic Raman 
transitions [17] or spin-resonance techniques [IB]. In this paper we focus on optical 
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transitions between one of the ground states, |1), and some excited state, |e), which 
may be short-lived. We define the optical raising and lowering operators as 

<7+ = |eXl|, <r- = |lXe| . (3) 

We presume that our light is completely ineffective at inducing transitions between |0) 
and |e) either because |e) is too far off-resonance, because of a prohibitive selection 
rule, or some combination of the two. One example of such a system is provided by a 
semiconductor donor-bound impurity, where the qubit states are provided by electron 
Zeeman sublevels and the excited state is provided by the lowest bound-exciton state. 
Other examples include the hyperfine structure of trapped ions. In this paper we will 
always refer to the matter qubit as an atom although it may be a semiconductor impurity 
or quantum dot comprised of many atoms. 

In an ideal case, off-resonant light results in the effective Hamiltonian 

H = -JnZ/2, (4) 

where n is the photon number operator. (An arbitrary state-independent optical phase 
term has been removed from this Hamiltonian). This Hamiltonian is not physically 
exact, but rather a desired result of a strictly dispersive atom-cavity interaction. Such 
a Hamiltonian is often assumed after the excited state |e) is "adiabatically eliminated." 

To motivate the next section, in which a more detailed model appears, let us first 
explain our entanglement protocol as if we had a lossless channel. Qubit 1 would initially 
be rotated into state (|0) 1 + {1)^/^/2. When a pulse of coherent light, which we call the 
"bus," reflects off of the cavity containing the qubit, the interaction lasting time Q\j J\ 
results in the unitary evolution operator 

U l = exp(i9 1 nZ 1 /2). (5) 

After this interaction, the bus with initial coherent state amplitude a and the rotated 
qubit would be described by the state 

^[|0>i + l«> = ^tlOKK^ 2 ) + |l»e-^ 2 >]. (6) 

After traveling the fictitious lossless channel, the bus would then undergo an identical 
interaction with a second qubit, which had been synchronously rotated into the same 
initial state as the first. The resulting selective phase shift by 9 2 = 9 1 = 9 yields the 

state 

' ' '0> 2 |a*> + (|0> 1 |1> 2 + |1> 1 |0> 2 )|«> + ll^ll^la-*)]. (7) 



2' 

If a is made infinitely large, the states \a), \ae td ), and \a~ ld ) are orthogonal, so that they 
may be distinguished by some detection scheme such as homodyne detection. Employing 
such a scheme and keeping only those events in which the bus carries zero phase shift 
(state | a)), we post-select the qubits into the maximally entangled Bell-state 

* + > = ^ f [|0} 1 |l) 2 + |l} 1 |0} 2 ]. (8) 

This Bell-state could now be used for quantum teleportation. 
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Unfortunately, the situation is more complicated in the presence of channel loss. 
In this case noise is introduced during the transmission of the bus, requiring a density 
operator approach to describe the state. In order to keep the noise reasonable, a finite 
value of a must be taken. Then the phase shift of the light cannot be perfectly resolved, 
and so a more detailed look at the homodyne detection is needed. We pursue this 
analysis in the following sections. 



3.2. Effective Interaction with Channel Loss 

In order to generalize and allow for the introduction of noise, we write the state of qubit 1 
as a general density matrix p 1 , possibly already entangled to other qubits. Equation |H1 
may then be generalized to 

Ux \a) pi (a\ U\ = IfciZ^p^ZJl, (9) 

where 

(3 X (Z X ) = ae iZl9l/2 . (10) 

To clarify the notation, Z\ is the operator of equation (j2J) operating on p\\ the meaning 
of operators inside kets is unambiguous in the basis where these operators are diagonal. 

In the state described by equation (jOJ), the qubit may be highly entangled with the 
bus. The maximum degree of entanglement can be found from the entropy increase after 
tracing over the optical states; as a — > oo, this approaches 1 — exp(— 4\a\ 2 sin 2 #i/2) bits, 
indicating that the entanglement increases as we increase the average photon number 
of the pulse \a\ 2 . Unfortunately, in a real system, as a is increased so does the amount 
of quantum information leaked to the environment by various forms of loss. In order to 
estimate the amount of entanglement actually available from this interaction, we must 
address these non-idealities. 

During the dispersive atom-cavity interaction, there is some probability that the 
atom is brought to state |e), after which it may undergo spontaneous emission or undergo 
non-radiative decay. This probability leads to a very damaging decoherence, since 
a single spontaneously emitted photon, phonon, or Auger-ionized particle reveals the 
"which-path" information of the qubit. In the idealized picture governed by equation (J3J), 
such a form of loss is not present, but we will consider it in detail in ensuing sections, 
where we refer to such processes as internal loss. 

For long-distance entanglement distribution, the dominant source of loss of photons 
is not from the cavity but rather from the ensuing communication channel. This loss 
depends on the length of the fiber connecting the cavities, although it may also be 
affected by imperfect mode-coupling to the cavity and the need to shift the wavelength 
of the cavity emission to a more convenient telecommunication wavelength PHJEO]- We 
refer to such processes as external loss. 

For now, let us consider only the dominant source of loss: external loss from the 
fiber and associated mode coupling. In this case the bus coherent state is damped from 
\f3i(Z 1 )) to \/T f3\{Zi)\ where T is the total transmission from one cavity to the next. 
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The lost photons are in principle lost into many modes. Although the physics of each 
loss mode may be complicated, we may treat it as we would a series of beam splitters. 
We write the lost photons in state 

\L) = Y.D m {a ( t ) e lZiei/2 ) |vacuum) , (11) 

m 

where D m (affl) is the displacement operator for loss mode m, and the total power 
in these modes must sum to the total amount of lost optical power, I] m |a:^| 2 = 
(1 — T)|a| 2 . We may immediately trace over these modes. As a result of a single atom- 
cavity interaction followed by loss, the effective interaction accomplishes the quantum 
operation 

\a) Pl (a\ -> \VTf3 1 (Z 1 )}Q 1 (p 1 )(VT(3 1 (Z 1 )\. (12) 

The superoperator Q\ is given by a rotation and a phase flip with probability A_: 

Qi(p) = e iZ ^ 2 [X +P + X^Z.pZ^e-^ 2 . (13) 

The parameters A± are given by (l±e _7l )/2. The losses 71 and phase shifts £1 are given 
by 

7l = |a| 2 (l-T)(l-cos^), (14) 
& = |a| 2 (l-r)sin0i. (15) 

The rotation by angle £1 can be quite large for realistic transmission T, but such a 
single qubit operation may be undone locally. A key assumption here is that this phase 
is accurately known, which will require prior knowledge and stability of the angle 9\ 
and real-time measurement of a and T. The loss of coherence going as exp(— 71) is the 
dominant and unavoidable source of the reduction of entanglement fidelity. 



3.3. Long-Distance Entanglement Distribution 

We now show how to use this semi-ideal interaction to distribute entanglement over 
long distances, referring to Fig. El The goal is to post-select a bus state that lacks 
"which-path" phase information from its two interactions with two qubits. 

Unfortunately, the two qubits are likely to have slightly different interaction 
constants J. The qubits may be somewhat inhomogeneous, but even perfectly 
homogeneous qubits will give different angles. As we discuss in section the angle 
depends on a, which is reduced due to fiber loss, and on the pulse length, which is 
increased due to fiber dispersion. Such an angle difference may be compensated as 
follows. After reflection from the first cavity, the bus pulse and a local oscillator (lo) 
pulse are transmitted nearly simultaneously down a single fiber. Note that the LO pulse 
undergoes the same random phase shifts that the bus pulse may have accrued as it 
traveled the fiber. The LO pulse is then immediately split in order to accomplish a 
small displacement of the bus pulse prior to interacting with the second cavity. Such a 
displacement occurs by mixing the very strong LO pulse, with amplitude a L0 , with the 
bus pulse at a 1 — x : x beam-splitter, where 1 <C 1. One output of the beam-splitter 
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Figure 3. Schematic for long-distance entanglement distribution. A laser pulse is 
split into a local oscillator (lo) pulse and a probe pulse. The latter reflects from the 
left cavity, travels nearly concurrently with the LO pulse, and reflects from the right 
cavity. Homodyne detection is then performed to post-select the entanglement. 



gives the bus displaced by a term of order xa L0 , while the decoherence caused by the 
light dumped at the other output is of order x 2 a LO , a term which adds to the external 
losses already discussed. Now, the displaced bus pulse interacts with the second cavity 
and qubit, this time resulting in phase ±6*2 for each eigenvalue of the Z 2 operator acting 
on the state of the second qubit, p 2 . This displacement and interaction yield the state 

U 2 D((3 T )\VT(3 1 (Z 1 )}Q 1 (p 1 ) <g> p 2 (Vf(3 1 (Z 1 )\D\(3 T )U 1 2 = 
e iMZi)\ R + i^Q^p,) ®p 2 {R + il\e~ ich ^\ 



(16) 



where 

R + il = (3 T e iZ ^ 2 + VTae^ +z ^ 2 , 
<t> T (Zx) = lm{VTae- iZl/d2 (3 T }. 

The ideal choice for the "tuning-displacement" amplitude /3 T 
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With this choice, we find 
R(Z 1 ,Z 2 ) = 

I(Z U Z 2 ) = 



sin 9 2 /2 



I — $1 02 ( 9 $2 

via sin — sin — cot — 
2 2 V 2 

VTa sin — cos + ^2 



Z\Z' t 



;i7) 
:is) 



(19) 

(20) 
(21) 



The term exp[i(f) T (Z 1 )] represents a Z-rotation of qubit 1 by angle —\v / Ta\ 2 sm(9i/2 — 
6*2/2) sin(0i/2)/ sin(0 2 /2). Like the rotation by £1 discussed in the last section, this large 
single qubit rotation may be removed, and so we will consider it no further. 

To generate entanglement, we note that / = in the subspace where Z\ + Z 2 has 
eigenvalue 0. If the initial state of the qubits is (\0) 1 + |l) x ) ® (|0) 2 + |l) 2 )/2, which is 
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achieved by simultaneous, phase-coherent 7r/2 Y-rotations of the two distant qubits, and 
if we project onto the Z\ + Z 2 = subspace, we achieve the maximally entangled Bell 
state \^ + )- Such a projection can be achieved by measuring a phase of zero from the 
bus state. In a probabilistic picture, this phase shift of zero implies one and only one of 
the atoms interacted with the pulse, but without information as to which atom caused 
the interaction, the quantum state is described by a superposition of both possibilities. 

Of course, the phase of the pulse is a continuous quantum variable, and it may 
only be weakly post-selected with finite probability. Such weak post-selection may be 
accomplished by p-homodyne detection, in which we interfere the bus state with a LO 
pulse 7r/2 out of phase from the initial bus pulse %. The difference photon number in 
the two output ports is measured; the result is proportional to the projected eigenvalue 
of p = (a — a^)/2i §. As a result of this measurement, the qubits are projected into state 

(p\R + U)Qx( Pl ) ® P2 (R + tI\p} = U(p)G(p)Q 1 ( Pl ) ® P2 G(p)U\p), (22) 

where 

U(p) =exv[-i(2p + I)R], (23) 

G(p) = (2/7r) 1 / 4 exp[-(p + /) 2 ]. (24) 

How close is this conditional state to the desired Bell state |\E' + )? Since (Z\ + 
Z 2 )\^ + ) = and Z 1 Z 2 \^ + ) = - |^ + ), we see that |^ + ) is an eigenstate of U(j>), 
and so U(p) has no effect on the fidelity with respect to |^/ + ). We therefore focus 
attention on the Gaussian projection operator G(p). The projection described by this 
operator may be understood graphically in figure 0] From this picture we see that any 
measurement with result near p « projects the density operator into the parity-odd 
Zi + Z 2 = subspace. To post-select these cases, we keep only measurement results 
where — p c < p < p c . The probability of such an event is 

P s = f Pc dp Ti{G 2 (p)Q 1 ( Pl )®p 2 } 

(25) 

2 erf [y/2p c ] + erf [y/2(p c + y/Td)} + erf [y/2(p c - VTd)} V ' 

~ 4 ' 

Here we have used the important parameter 

d = 2a sin(0i/2) cos(# 2 /2), (26) 

which we refer to as the distinguishability. Note that if Q\ = 9 2 = 6, then d = a sin 9. 
The fidelity of the post-selected state varies depending on the measurement result p. 
If we discard the information of the precise value of p and automatically keep only 

| Other measurement schemes such as x-homodyne [2] or photon-counting [22 are possible, but 
for high fidelity in the presence of large external losses and high probability of success with realistic 
detectors, p-homodyne shows improved performance. 

§ This definition of p is convenient for the present calculations, but is not the only convention; 
in particular this convention leads to the commutator [x,p] = i/2. Previous presentations of 
hybrid architectures for quantum information processing, such as |21l I22j , have used the convention 
p = i(at — a), in which case [x,p] — 2i and the definitions of U(p) and K(p) are altered. 
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Figure 4. A quasi-probability distribution function (e.g. Wigner function) is shown 
on the left. The p-homodyne measurement yields a probability distribution function 
(G 2 (p)), shown on the right, found by integrating the Wigner function over x. The 
expectation value is evaluated for different choices of qubit states. The Gaussian peaks 
corresponding to finding the two qubits in states |00) or |11) are displaced from the 
p = origin by ±^/Td due to the state dependent phase-shifts totalling to ±9 that 
accrue during the interaction. However, states 1 01) and 1 10} accrue no total phase shift 
at the end, so if p is measured between —p c and p c , this subspace is approximately 
post-selected. 



those instances where p falls in the post-selection window \p\ < p c , the resulting density 
operator is calculated as an average of possibilities, 

A = ^ r d P U (P) G (P)QM ® P2G(p)U\p). (27) 

The fidelity with respect to |^/ + ) is then 

^ + ) (28) 

(l + e -7i) erf [72p c ] 

" 2 erf [V2p c ] + erf [y/2(p c + VTd)} + erf [y/2(p c - \fTd)] ' 

If c£ — > oo, then the postselection operation G(p) works very well, and p± 2 is entirely 
contained in the desired Z\ + Zi = subspace. However, the fidelity reduction due to 
fiber loss is characterized by 71 = d 2 (l — T)/2 in the low 9 limit. As d — > 00, then, the Q\ 
super-operator casts the two-qubit density operator into a classical superposition of |01) 
and 1 10), rather than the desired |^/ + ) state. For quantum coherence to be preserved, d 
must not be made too large. 

Figure El shows the fidelity of the final two-qubit entanglement as a function of 
the distinguishability d and the post-selection window for T = 0.67. We see that 
at each p c , there is a maximum fidelity at which d is large enough to allow good 
post-selection but not so large that fiber loss destroys all coherence. This maximum 



P12 
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0.00-1 , . , , 1 

1 2 3 4 5 
d = a sin 9 

Figure 5. The top figure shows the fidelity with respect to Bell state | V I /+ ) as 
a function of distinguishability d for linearly varying values of the post-selection 
measurement window parameter p c between and 2.0. The bottom figure shows the 
corresponding probability of success. The calculation assumes T = 0.67, corresponding 
to 10 km of telecom fiber and ^ = ^ « !■ 




Figure 6. (a) Initial entanglement fidelity F = ( l J /+ |/9i2| v I' + ) vs. probability of 
success. Each curve is labeled by the distance between qubits, £, normalized by the 
attenuation length of the fiber, £q. The probability of success is increased by increasing 
the post-selection window parameter p c , and d is chosen for each p c to maximize the 
fidelity, (b) The values of the distinguishability d and post-selection window parameter 
p c which lead to the maximum fidelities of subfigure (a). 
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fidelity is largest at p c = 0, but this condition means that the scheme would succeed an 
infinitesimal number of times. At larger values of p c the probability of success increases, 
but since states outside of the desired state become more probable the optimal 

fidelity decreases somewhat. This optimal fidelity as a function of the probability of 
success is shown in figure a), for several different possibilities of fiber length. The 
range of distinguishability d required to achieve these curves is shown in figure EJb) • 
The calculations leading to these curves assume that 9 <C 1, in which case the results 
depend to a good approximation only on d and not independently on a or 9. 

In figure |3 the fiber length between repeater stations, £, is measured in units of the 
attenuation length £ [i.e. T = exp(—£/£ )]. For fused-silica fibers at telecommunication 
wavelengths, the attenuation length £ can be as high as 25 km. It is clear that much 
higher fidelities are available for communication distances much shorter than this length, 
but then the spacing of repeater stations may be too small to be practical. For longer 
distances, the reduced fidelity will require more purification ||. Optimizing the distance 
between repeater stations will involve a large number of trade-offs, and the best choice 
will depend on the efficiency of the purification protocol. 

As an example, a typical working condition is £/£q ~ 0.4, corresponding to 10 km 
of fiber and T=0.67; these choices were used for figure 03 If we choose p c = 0.5, 
the probability of success is 36% and the initial fidelity is 77%. Since the light 
source can be a normal stabilized laser and detection is extremely efficient, this 36% 
probability of success is not degraded by source or detector efficiency. The rate of initial 
entanglement generation is therefore extremely fast in comparison to most other schemes 
for entanglement distribution. 

3.4- The Measurement- Free C-Z Gate 

The rate of final, long distance entanglement will depend on the efficiency of the protocol 
for entanglement purification and swapping , which depends on the fidelity of local 
operations (especially two-qubit gates such as controlled-NOT). A complete architecture 
for a quantum repeater requires some way of achieving these local operations. 

The entanglement distribution scheme discussed in the previous section could in 
principle be used for these local operations, but its probabilistic nature leads to very 
inefficient purification protocols. Deterministic local operations not involving post- 
selection are possible using the same resources we have already assumed for entanglement 

|| Longer station-to-station distances are possible if the detectors are placed at a mid-point equidistant 
from the two stations. In this scheme, four pulses are used; a bus and LO pulse are sent independently 
from each of the two qubit stations to the central detector station, which corrects for phase differences 
in the two paths, interferes the two bus pulses at a beam splitter, and performs a pair of homodyne 
measurements. One homodyne measurement makes a p-projection for entangling the qubits, while 
the other makes an x-projection to disentangle the output of the second port of the beam splitter. 
With this geometry, a similar analysis yields that the parameter 71 is doubled while the fiber length 
I is halved. For a total qubit-to-qubit distance of 10 km and 36% success probability, as used in our 
previous example, the fidelity only increases slightly to 78%. However, in this case the fidelity remains 
above 50% at 50 km with up to 20% success probability. 
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distribution, i.e. the effective interaction of equation (j3J) with d ~ 1. Ways to achieve 
quantum logic gates with these resources were discussed in One method uses this 
effective interaction to approximate controlled displacements, and employs the Berry- 
phase accumulation under those displacements to achieve a controlled-sign (G-Z) gate. 
The specific choice of operations to construct the gate is as follows. The gate occurs 
between two qubits, each in separate cavities connected by a short, local waveguide 
network. First, a dispersive interaction with qubit 1 is performed, as described in 
section 13.21 Then the optical bus in state is mixed with a LO pulse at a 

beamsplitter with amplitude and phase chosen to achieve a displacement of a(i — 1). 
Following this displacement, an effective interaction between the bus and qubit 2 is 
performed. The bus state may now be described as \/3 2 (Zi, Z 2 )) , where 

P 2 (Z U Z 2 ) = e iZ ^l\p x {Z x ) + (z - IK]. (29) 

We have written 8 2 for this second interaction in case the controlled rotation during this 
interaction is slightly different from the first. After this interaction, the bus is again 
displaced, this time by — a(l + z), after which it interacts with the first qubit again. 
This brings the optical bus to state |/3 3 ), where 

P 3 (Z U Z 2 ) = e iZ ^ 2 [j3 2 (Z li Z 2 ) - (z + l)a 2 }. (30) 

The cycle is completed by a displacement by a(l — z) and a final interaction with the 
second qubit, resulting in state |/? 4 ) , where 

P,(Z U Z 2 ) = e iZ ^ 2 [fo{Z u Z 2 ) - (i - l)a 3 ]. (31) 

If 9 n is the same for every interaction, then this final coherent state amplitude is 
approximately a constant term, — i\a\ 2 , plus a term proportional to d 2 involving 
the product operator Z\Z 2 . (If 6 n differs during each interaction, the gate requires 
modifications of the displacements to compensate, similar to the "tuning-displacement" 
used in the previous section.) The C-Z gate results from the phase accrued when tracing 
over the bus state, in addition to the phases that accrued during the displacement. The 
final state-dependent phase contains terms proportional to Z\ and Z 2 , which correspond 
to single qubit rotations, but also the term \a8\ 2 ZiZ 2 /2, for small 9. With the bus state 
removed and single qubit rotations corrected, this remnant phase is a nonlocal unitary 
operator acting on the two qubits, which may be used as a C-Z gate if \a6\ 2 = ir/2. 
More details of this calculation, including the effects of optical loss, may be found in 
Appendix A 

As in the scheme for entanglement distribution, the performance of this C-Z gate 
is limited by both internal cavity losses and external losses. As this is a local gate, 
these external losses are dominated by the interfaces between short waveguides and the 
required cavities and beamsplitters. 

The total fidelity of each gate as a function of local, external loss is shown in figure [7| 
This fidelity reduction should be compared to the error due to internal loss. The C-Z 
gate we analyzed interacts with each of the two qubits twice, and therefore the total 
error due to internal losses goes as the calculated fidelity of the atom-cavity interaction 
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Figure 7. This shows 1 — F, where F — Ao is the fidelity of the measurement-free 
C-Z gate, as a function of the local external optical loss, for several different values of 
9, maintaining the requirement ct6 ~ 1. See |Appendix A] 



to the fourth power. Whether internal or external losses will dominate for this gate will 
depend heavily on the details of the system used to implement it. 

3.5. Final Communication Rate 

To estimate the final communication rate for this proposal, we revisit the entanglement 
purification and swapping protocol discussed in section El We repeat the same 
simulation, except using more realistic noise models. We use the calculated initial 
density matrix given by equation (|27j) for the initial state, estimating that external 
loss is dominated by the attenuation length of standard fused silica fiber at telecom 
wavelength. We use the noise model developed in |Appendix A to model gate errors 



during entanglement purification and entanglement swapping; this noise is dominated 
by local loss in the short waveguides between qubits in a single repeater station. We 
assume perfect single-spin rotations. As an illustratitve example, we again presume 
N + 1 = 129 repeater stations, each spaced by 10 km, with 2 + 2 log 2 A^ = 16 qubits in 
each station. Each qubit is in its own cavity. Such an architecture could be implemented 
using fiber-optic waveguides and switches, or on a single chip using planar photonic 
crystal waveguides and switches. 

A typical number of purification steps under these conditions might be three rounds 
of entanglement purification of nearest-neighbor, initially generated entanglement, two 
rounds of entanglement purification after the first few levels of entanglement swapping, 
and one or zero steps of entanglement swapping for the last few levels of highest- 
distance entanglement swapping. The exact number of steps is varied from simulation to 
simulation to achieve different communication rates and different levels of final fidelity. 
Figure |H1 shows the final communication rate and fidelity for several example choices 
as a function of local external optical loss. Communication rates approaching 100 Hz 
appear to be possible. 
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Figure 8. Rates of final generation of distant (1280 km) entangled pairs resulting 
from Monte-Carlo simulations of the nested entanglement protocol with 16 qubits per 
station and 127 intermediate stations separated by 10 km. The simulations run until 
5 pairs are generated; the average time and standard deviation between generated 
pairs is plotted. The error modes included in the simulation are the initial fidelity 
reduction due to external loss in the 10 km fiber and distortion due to local external 
losses during the measurement-free C-Z gate, as described in |Appendix A| Different 
numbers of purification steps are used to achieve different final fidelities. 



4. Calculation of the Dispersive Light-Matter Interaction: Methods 

Up until now, we have assumed the effective Hamiltonian of equation (@J), and we 
have argued that this is a sufficient interaction to achieve a fast, long distance 
quantum communication protocol. Unfortunately, the idealized interaction described by 
equation (j3J) is not physically available in a CQED system, except as an approximation. 
This approximation assumes the "bounce" fidelity of the dispersive interaction is 
negligibly close to one. In this section, we analyze the physical parameters of the 
atom-cavity system required to satisfy this assumption. 

Most existing theoretical analyses of CQED designs for entanglement generation 
are only appropriate in regimes different from those considered here. For the strongly 
coupled regime, evaluation of dynamics in a dressed-state picture is an effective means 
of calculation, but in the weak or intermediate coupling regime, the large widths of 
the dressed states make such calculation techniques inappropriate. Weak-coupling or 
intermediate-coupling calculations usually make one of several approximations. One 
simplifying assumption is that only single photons or weak-coherent states interact with 
the cavity, heavily limiting the dimensionality of the equations to be solved. This has 
been effective for studies of spontaneous emission and the proposals for CQED devices 
in the weak-excitation limit, but it will not be effective here. Typically, when more 
photons are introduced into the cavity, a numeric approach is required. For very 
large photon numbers, a full-quantum analysis can be computationally intensive; an 
appropriate approximation is the semi-classical optical Bloch equation approach. We 
will primarily work under such assumptions in the present study. 

We establish notation by first describing and solving the empty cavity problem. We 
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then compare the interaction in a loaded cavity using an interaction picture that only 
considers dynamics in a frame with the empty-cavity dynamics removed. We begin with 
an analytic, perturbative approach, which expands the dynamics in the atom-cavity 
coupling g and the bare-atom relaxation rate r _1 . These are presumed to be much 
smaller than the time for light to leak out of the cavity, ujq/Q (weak or intermediate 
coupling regime.) It is not assumed that the atom is never excited, or that the photon 
number is small, although we will see that the expansion fails if the photon number 
is large enough to saturate the interaction. For more accurate results, we will employ 
a semi-classical simulation of the full many-photon atom-cavity dynamics. We will 
check the results of this semi-classical approach against the results of the perturbative 
calculation in the low photon number regime and against a fully quantum analysis in 
the high photon number regime. 



4-1. Empty Cavity 

The empty-cavity Hamiltonian is 

H = u ala + u\b\b x + « K ^( b x a o ~ b x a o)- (32) 

A A 

We are working in a frame rotating at the frequency of the |e) <-> |1) transition, iv a . 
The operators 6a (b\) annihilate (create) photons in the eigenstates of the waveguide 
to the cavity, with eigenenergies h(cu a + The cavity- waveguide coupling constants 
K\ are assumed to be real. The sum includes lossy modes and absorptive losses in the 
cavity mirrors. The operators ao and a\ respectively annihilate and create photons in 
the single cavity mode. We further make a transformation to wave-packet modes 

« m = E & K m , (33) 

A 

for m > 0, where U\ m is a unitary matrix. (When summing over modes in the following 
formalism, the m = cavity mode is also included.) The element u\ rn is the Fourier 
component of a traveling pulse with label m. Two modes are of particular interest: 
a m annihilates a photon in the wavepacket incident on the cavity, and a 0VT annihilates 
a photon in the wavepacket of the output mode which couples to the waveguide for 
communication. Hence u\ tlN describes the input pulse shape f m (r,t) according to 

u x , m = / d 3 f / dt ^(r^Ur, t), (34) 

JV J— 00 

where ip\(r) is the spatial shape of the wavepacket for waveguide eigenmode A. 

The solution of the empty cavity problem is well known; here we treat it as a 
scattering matrix for the operators a m . The Heisenberg equations of motion are 

a (t) = -iu a (t) + K X b\(t), (35a) 

A 

b\(t) = -iu\b\(t) - K\ao(t). (356) 
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Using Laplace transforms, a(s) = e st a(t), we arrive at the scattering matrix 
equation 

h(s) = J2 S ^m( S ) a m, (36) 

m 

where 

SooO) = , ■ 1 , jk, (37a) 

s + itu + 7/2 

Sxo(s) = -S 00 (s)^—, (376) 

Som(s) = -J2S X0 {s)Uxm, (37c) 
A 

c / v U\m — K\S 0m (s) fnnj\ 

b\ m {s) = — . (37 a) 

S + lUJx 

We assumed a sufficiently broadband spectrum of output coupling and absorption modes 
to lead to an exponential cavity decay function [23J; this came from the mathematical 
association 

H^- = l (38) 

implying that any optical power in the cavity leaks out of the cavity as e _7i . 

We also introduce an average output coupling factor re (without subscript) as 
follows. Suppose at t — the cavity contains a sing le photon, |^(0)) = 4|0). Then 
after some time the photon leaks into all possible modes A, 

lim|tKt)> = lim£S o t A (t)&t|0). (39) 

A 

Many of the modes indexed by A will not be included in the desired output wavepacket. 
We therefore seek the overlap of this state with that of a single photon in the input 
mode, \4> w (t)) = J2\ e ~ luJxt b\u\ tlN |0) in the limit where t — ► oo. We arrive at the overlap 
integral 

m) I M*)) = E v KXU T To = ^in(-^o - 7/2). (40) 
x H^a - ujo) - 7/2 

where F w (s) is the Laplace transform of the input pulse shape for mode m as it couples 
into the cavity. Two effects are present in this overlap: first there is the overlap of the 
pulse transform with the cavity filter function, indicated by the Laplace transform of 
the input light F m (s) at s = —iuj$ — 7/2. Second, there is the output coupling factor k; 
only a fraction of roughly re/7 of the light makes it from the cavity to the output mode. 
Using these definitions, we write 

SoAs) = J U f /9 - (41) 

S + IUJq + 7/2 

This function, which represents the convolution of the input pulse with the filter function 
of the cavity, will be used heavily in what follows. The pulse reflected from the cavity 
can be seen from equations PTj) and ()37dj) to have components 

S\, m (s) = — — — u x , m . KX j-<jKF m (s) , (42) 

s + lUJx V s + iu + 7/2 J 
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and therefore 

?<a.»w = 1tC-jti^ |f ™ w|2 (43) 

The first term describes the traveling wave corresponding to mode IN, the second term 
describes the interference from light that coupled in and then back out of the cavity. If 
the cavity mirror passed all modes with equal coupling, then the light entering the cavity 
would have the same shape as the input pulse, i.e. |F IN (p)| 2 = J2x |^a,in| 2 /(p + In 
this limit, we would have 

- r s + iuo + ^/2-K] . 44 > 



X!' U in,A'S'a,in(' s ) — 



IN 



x S + IUJ\ 



s + iujQ + 7/2 



If we further assume the cavity to be strongly overcoupled [k = 7), we arrive at the 
solution for reflection from an empty cavity found in [23J. 

The empty cavity is a linear scattering problem and may be solved with this 
scattering matrix approach. However, the presence of an atom in the cavity transforms 
it to a nonlinear problem. For this we consider the characteristic function 

X (ri,t) = Tr{D OVT ( V ,t)p(t)}, (45) 

where D 0VT (r), t) is the displacement operator for the time- dependent output mode OUT: 

D vt(v, t) = exp[r/4 UT (t) - r]*a OVT (t)}, (46) 

for cto UT (t) = J2x b\S Xm (t). If the input wave-packet is a coherent-state with amplitude a 
and the cavity is empty, we find 

X„M) = e-'^IlTr {e^^.W^ |5 A , IN (t)a)(5 A , IN (t)a| A e^,N(^a . (47 ) 

A ' 

Now, we may use the unitarity of the scattering matrix, i.e. 

£^,a(*)<V(*)=<Wn, (48) 

A 

to see that 

Xo ( v ,t) = e-l^l 2 / 2 ^*-"*", (49) 
the characteristic function of a coherent state in traveling wavepacket-mode OUT. 

4-2. Atom- Cavity Interactions in the Interaction Picture 

If the cavity is not empty, the time-dependence of p will be more complicated. In this 
case, we consider the density operator in the interaction picture, 

X (v,t) = Ti{D 0l]T ( V ,t)p(t)}, (50) 

where 

p(t) = e mot p{t)e~ iHot (51) 
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and 



ew(- l -^f + J2vb[(t)S x ^t) - h.c. 



\m 

2 

I 12 x (52) 

e -W 2 /2+ml-v*a m = jDiN ( r?;0 ). 



"H.c." refers to Hermitian conjugate. 

To analyze the effect of an atom in the cavity (which we call an atom, although it 
may be a quantum dot or impurity complex), we consider unitary dynamics governed 
by a Jaynes-Cummings term [23] . 

Wint = 9(%o- + + ajcr~), (53) 

as well as non-unitary atomic relaxation processes. In the interaction picture, this 
Hamiltonian becomes time-dependent: 

Wtat(t) = 9j2[ S o m (t)a m a + + alSl (t)o--}. (54) 

m 

We assume that relaxation is limited by lifetime effects, so that the full Liouville-von 
Neumann equation may be written 

f t =^[n- mt (t),p(t)] + mt)]- (55) 

The latter terms represent atomic relaxation at zero temperature. Considering a master 
equation approach, we use 

= -^-W + °-p(t) + ~pit)a + a- - 2a~ p(t)a + }. (56) 
Zt 

Here, r is the lifetime of the atom in the absence of the cavity, including both 
spontaneous emission and non-radiative decay, such as Auger decay of bound excitons 
in silicon. In principle, pure decoherence terms could be added to the above; however 
in the present study we neglect decoherence between the two ground states and assume 
the optically created coherences are strictly lifetime broadened. 

As we will see, the most important parameter for quantifying the performance of a 
particular atom-cavity system is related to the Purcell factor 

where r r is the spontaneous emission lifetime outside of the cavity. We will be considering 
systems where the atom is far detuned from the cavity frequency, so that the Purcell 
effect is weak. However, the key parameter quantifying the suitability of an atom-cavity 
system for strong dispersive interactions with low absorption is 
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This parameter is sometimes known as the "cooperativity parameter" or the inverse 
of the "critical atom number." We will see that if a, the center frequency of the 
pulse, and the pulse width are optimally chosen, the final fidelity of the dispersive 
interaction approximately decreases as exp(— d 2 /$), so a large $ is critical for high- 
fidelity operation. 

We now have enough formalism to clarify a terminology we have already employed. 
We say that a cavity is in the strong coupling regime if g > 7 and g > 1/r. The strong 
coupling regime is not optimal for this quantum repeater architecture. The optimal 
regime is the intermediate coupling regime; specifically this means that g < 7 but 
$ > 1. 

We estimate the atom-cavity coupling from the formula 

2 _ 3 u; a A 3 

where V is the mode volume, r r is the strictly radiative lifetime, A and u) a are the 
wavelength and frequency of the atomic resonance, and n is the index of refraction of 
the host material. This formula assumes an optimally aligned dipole at the antinode 
of the cavity field. We even apply this formula to silicon, where the indirect, phonon- 
mediated transitions are not simple electric dipole transitions. The slightly different 
probabilities for phonon absorption and emission means that the silicon system is 
not exactly described by the Jaynes-Cummings model, but these corrections are not 
expected to be important for modeling far off-resonant dynamics. Using equation (JoTJj) . 
then, and assuming a pulse resonant with the cavity, we have 

$ = JL x 1 x 1 x ^- x Q. (60) 
Att z 7 r r n a V 

The four terms of this expression will serve to summarize the importance of each physical 
parameter. The factor indicates the degree that light leaking from the cavity leaks 
into the desired output mode. For the present quantum repeater application, it is best 
that the cavity is overcoupled; i.e. that cavity loss is dominated by transmission through 
the output mirror, so that k/j — > 1. The factor r/r r is the ratio of the total lifetime 
of the atom divided by the radiative lifetime, i.e. the internal quantum efficiency of the 
emitter. This factor is detrimental in the case of donor-bound excitons in silicon, for 
example, where r is dominated by Auger recombination. The factor \ 3 /n 3 V indicates 
the importance of a microcavity in achieving a high value of $. Finally, the total cavity 
Q = uj & I 7 is a critical parameter; it must be large enough to compensate for deficiencies 
in the other factors. 



4-3. Analytic Approximation of Unsaturated Phase Shift 

We first present the phase shifts and internal loss parameters due to the atom-cavity 
interaction expected from an analytic, perturbative approach. This approach only 
provides limited value in the regime of high coherent state amplitude a, but will help 



Hybrid Quantum Repeater Based on Dispersive CQED 



24 



motivate the numerical calculations which follow. The details of this approach are 
presented in Appendix B 



To second order, we find that after the bus pulse completes its reflection from the 
cavity loaded with an atom in state |1), the coherent state remains a coherent state 
with amplitude a(l — i6 2 — L 2 ) . (Here, the subscript refers to the order in perturbation 
theory. Also note that this description, in which the atom in state |0) results in no phase 
shift and an atom in state |1) results in phase shift 9, differs from the description given 
by equation (jlj) by an overall optical phase shift.) The second-order state-dependent 
optical phase-shift 9 2 is 

a 2 f duJ \F m (-tu)\ 2 

9 2 = g k p.v. / - 61 

J 27T a; [(a; — u ) 2 + 7 2 /4J 

The coherent state amplitude is also reduced in this order with internal loss parameter 

g 2 ^ \F a 



For a narrow-band pulse which is off-resonant with the atom, we approximate 
\Fm(—iuj)\ 2 « 8{u — Up), in which case we find that a sees a small-angle phase shift of 
magnitude 

_ Q Ki 1 Ki , r \ 

02 ~ rj- 2~Ta\ = ^(Wr - U ). (63) 

Up[{Up -UqY + 774J UpTr'J 

Here we see the first simple principle which will be important in the design of 
systems employing this interaction. The largest phase shifts are available when u P = Uq, 
that is, when the center frequency of the pulse is on-resonance with the cavity. This 
makes sense; the more passes the light makes in the cavity, the stronger the dispersive 
interaction, and the light makes the most number of passes on resonance. For the 
remainder of this paper, we will assume this simple condition. 

In general, the phase shift might increase with smaller offset u P between the pulse 
and the atomic resonance. However, the smaller this offset, the larger will be the 
loss. This is already evident in the L 2 terms, but even if |F IN (0)| 2 is neglected, atomic 
dephasing will dominate the loss at small u P . This is seen in third order, under the 
off-resonant, narrow-band pulse assumption, in which we find 

— ) 02- (64) 

(A correction to this third-order loss occurs proportional to |F IN (0)| 2 and the first 
derivative of \F m (—iu)\ 2 at u — 0.) This shows a second important but simple principle: 
the lowest order loss term is an extra factor of the offset from the dispersion. This loss 
term is related to a reduction of the final fidelity of the interaction. 

If we go to higher orders in the perturbative expansion, we see a number of expected 
terms. In fourth order we find that — Q\ should be added to — i6 2 and similar higher 
order factors of the loss, corresponding to the expected behavior of a — > aexp(—i9 2 ). 
This suggests that our expansion is valid as long as 9 2 <C 1, which is indeed the regime 
of interest. However, already in fourth order, a term appears of the form — |ct| 2 #2, 
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suggesting that if a is made too large, a saturation effect occurs. As we will see later, this 
is definitely the case, and this means that the criterion for validity of the perturbative 
expansion is la) 2 ^ <C 1. Unfortunately, the very premise of the schemes discussed 
in section El requires that \a9 2 \ ~ 1. Hence, this perturbative approach is inevitably 
limited for further analysis, and accurate calculations of the fidelity will require another 
approach. 

Before leaving perturbation theory, however, we might note the order in which 
it predicts a few more important effects. The first term in which the coherent state 
amplitude is affected by atomic population in the excited state |e) is in fifth order. It is 
in this order that we must also begin to consider the lossy modes of the system, where 
the trace over other modes besides m =IN introduces new loss terms. It is not until 
sixth order that we start to see any non-Gaussian features in the characteristic function. 
This suggests that treating the light as a coherent state throughout the calculation is 
an excellent approximation, a suggestion that we will employ and then test numerically. 
Only at high values of a where the interaction is nearly saturated does a quantum 
treatment of the light become important. In this regime the fidelity decay due to 
internal loss is fairly low anyway, and therefore this regime should be avoided for the 
desired interaction. 

4-4- Optical Block Equations 

For a more accurate calculation of phase shifts in the presence of large values of a, 
we numerically solve the quantum master equation. Our approach is as follows. We 
first derive the master equation in a fully quantum setting in which any quantum state 
of light is allowed. In the regime of interest, we will find that the light may always 
be described as a coherent state. Assuming this condition, we enter the semiclassical 
approximation of the optical Bloch equations. 

4-4-1 ■ Quantum Master Equation. In our interaction-picture approach, we are 
comparing the output pulse to that expected from an empty cavity. For this we only 
model coherent dynamics in the IN mode and the cavity mode. This problem is nearly 
equivalent to an analysis of the coherent interaction of the single cavity mode with a 
time-dependent coupling given by gS m {t). However, the lossy modes are important for 
calculating the fidelity, and in developing single-mode optical Bloch equations we must 
incorporate the effects of such loss. 

These lossy modes are incorporated with a master equation approach, in which we 
assume the density operator may be written p ® |0)(0| L , where the L subscript refers to 
photons lost from the cavity due to leaky modes (including absorption at the mirrors). 
As the system evolves, photons enter these leaky modes, but unlike cavity photons they 
are immediately lost, resulting in the assumption that the L subspace is roughly vacuum 
at all times (Born approximation). Our master equation then results from 

Pit) = Tr L |-i[W(f),p(f) ® |0)(0|J - jf*df[W(f), [H(t'),p(t') ® |0)(0|J]} . (65) 
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The first term represents single-mode coherent dynamics. The second term is the term 
which introduces loss into leaky modes; it may be written as 

-g 2 J*dt'J(t,t')({o- + o--,p(t l )}-2o--po- + ^ +i)C{t,t')[a + a-,p(t% (66) 
where the real relaxation functions J and JC are given by 

j{t,t') + iic{t',t)= J2 s 0m (t)sUt f ) 

rn^lN (67) 

= (l-e-*)e- z ^-S ^t)Sl (t'). 

To simplify this master equation, we again take the limit where cavity and the pulse are 
both far off-resonant from the atom. In this case J~(t,t') + i)C(t',t) oscillates much 
more quickly than the time-dynamics of p. In this regime we may also make the 
Markov approximation. In the intermediate regime of interest here the Born-Markov 
approximation recovers the well-known Purcell effect. Only deep into the strong coupling 
regime does this assumption break down. 

As discussed in section EOl we presume that u P = u , i.e. that the pulse is resonant 
with the cavity (and both are offset from the atomic transition by u>o). Having clamped 
these two frequencies, it will now be convenient to work in a different rotating reference 
frame. In this case, instead of a frame rotating at the atomic frequency, we work 
in one rotating at the center frequency of the optical pulse (and the cavity mode). 
Correspondingly, we abbreviate 

S(t) = e^%, IN (t). (68) 

Then we may write 

J(t,t') + iJC(t',t) = e -^o(t-t')[ e -7(*-*')/2( 1 _ e - 7 t) _ s(t)S*(t')]. (69) 

In the Born-Markov approximation and in the long, off-resonant pulse limit we need 
consider the integral of this function over t' as t — > oo, for which we use 

g 2 J(t,t')^ -L^oMt-O, (70) 
ul + 7 2 K 

Since we assume the overlap of a spontaneously emitted photon with our pulse is 
negligible, the effect of the other modes to which the cavity couples is only to shorten 
the atomic lifetime. 

This single-mode approach is more readily tackled numerically with a c-number 
representation. In order to keep track of the atomic dynamics, we define several "partial" 
characteristic functions, defined by 

X ' k (v,t) = TT(j\D w (v)p(t)\k), (72) 

where states \j) and \k) are atomic states and the trace is over the optical field. With 
this notation, the complete system of c-number master equations for the characteristic 



g 2 JC(t',t) ->- ^_— 8{t-t'). (71) 
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function are 
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-2V X ee (v,t), (73a) 
+ 2T X ee (v,t), (73b) 
+ ( i n-T) X cl (v,t), (73c) 
(73 d) 
(73 e) 



Here 2r is the total decay rate of the atom in the cavity, including the influence of the 
Pur cell effect: 

1 + f(u P ) J_ 
• i 



2T 



(74) 



where r nr describes non-radiative decay % Also Q is the atomic detuning from the cavity, 
including the ac-Stark shift, 

1 



n 



uJ T . 



1 + — ^(Cjp 

7r r 



(75) 



Accurate numerical solutions of this system of equations at high values of a are 
computationally intensive, although we do so for some parameter sets as a check of the 
semiclassical approximation that forms the core of our results. For such calculations, 
we expand x^iv^t) m a truncated Hermite- Gaussian basis and solve the corresponding 
high-dimensional ordinary differential equation using Runge-Kutta techniques. For most 
of the parameter space discussed below, the semiclassical approximation is sufficient; 
the x^ k {v) calculated from a full quantum analysis are indistinguishable from the 
corresponding semiclassical results to the accuracy of the calculation. We only see 
appreciable non-Gaussian states of light at high values of a and low offsets when the 
atom-cavity interaction is highly saturated. Here a self-phase modulation effect occurs 
and the quantum phase-uncertainty becomes larger than expected for strictly coherent 
states. However, even in the semiclassical approximation, this regime is not appropriate 
for dispersive interactions with optimal fidelity, and so these quantum effects are of little 
consequence to the present study. 



4-4- 2. Semiclassical Approximation. The assumptions underpinning the semiclassical 
approximation are that the quantum state of the light during the interaction is always a 

If Our derivation of the Purcell effect leading to equation 174(1 omits several factors specific to the 
geometry of the cavity; for a more complete treatment, see |24j . for example. These factors have 
only small effect on the phase shifts in which we are ultimately interested, however, and so we use 
equation 1)74(1 in all simulations. 
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coherent state, and that it always remains unentangled from state |e). Then the density 
operator has the form 



(76) 



m = ia(t)X5(f)i ® [p c :(t)° + °- + ^(0) - P c :(t)]a-o-+ + P ?(t)o- + + ft®* 

+ \a(t))(*\®{pf(t)\e)(0\+pl°(t)\l}(0\ 

+ \a)(&(t)\®{pl c (t)\0)(e\+pl 1 (t)\0}(l\ 
+ |a)H®p°°|0)(0|. 

If we use this density operator in our equations for X \Vit) and focus on rj = 0, we 
arrive at the optical Bloch equations. The first three of these are 

pT = ig[S*(t)a*(t)p?(t) - S(t)&(t)pl c (t)} - 2Tpf(t), (77a) 
pf = igS(t)a(t)[2p™(t) - P l\0)} + (ifi - r)pf(t), (776) 

The only source of loss in these equations is from atomic decay. This may be seen by 
noting that these equations may be combined to derive, without further approximation, 

,2 ,„,2 Pl°(t) nl*f?V) 



Since p| e (t) — > as t — > oo, we have 



^^-wT^Mr (79) 

All optical loss is ultimately due to atomic decay. Optical loss from the cavity 
independent from atomic decay is already incorporated into the definition of S(t). 

From equations (77), we may approximately solve for the total phase shift and 
optical loss in the limit of a narrow-band, far detuned pulse. This approximation is 
obtained by first assuming gS(t)a(t) is constant in time, with value gSa, and solving 
the equations for p^(t) and p| 1 (t). (A similar approach was taken by [2]). These 
approximate solutions are 

oee(f] ^ (twrnai* 

Pa ( ^pTiiWS ( } " (81) 

Now, we presume this solution for p^(t) is maintained as S(t) varies in time. Then we 
integrate equation ()77c|) to find 

Q-iT 



a(t) ~ a 



ig 1 [ dt'\S(t') 
Jo 



(82) 



T 2 + tt 2 + 2g 2 \Sa\ 2 _ 

We emphasize that this approximation was not rigorously derived, and should be treated 
with caution. However, it shows the features expected from the previous section. 
First, in the t — > oo limit, and for the far-detuned, narrow-band pulse we have been 
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assuming, we estimate the integral over g 2 \S(t)\ 2 as $I\ Second, we assume that this 
approximation is the lowest order of an exponential solution to lead to our approximate 
equation 

lim -A2 « exp -z$r — ■ 09 ■ - .- . 83 

<-*oo a r 2 + n 2 + 2g 2 \Sa\ 2 J 

In the limit of vanishing a and T/Q = 0, this equation approaches 1 — i02, the result 
we obtained more rigorously in section 14.31 and |Appendix B| Further, if we look at the 
lowest order correction in T/Q, we obtain 1 — i9 2 — L$. This equation also correctly 
shows that as a is increased, a saturation effect occurs and the phase shift and optical 
loss both decrease, which begins to be evident in fourth order perturbation theory. The 
threshold a for which this occurs depends on the average value of S(t), which in turn 
depends critically on the pulse length. 

The remaining Bloch equations must be considered in order to calculate the fidelity 
of this operation, which is degraded by internal loss. These equations are 

pf = -iga(t)S(t)p™(t) - [i(A - O) + T + c(t)]pf(t), (84a) 
pl° = -igofS*{t)pf{t) - [iA + c(t)]pl\t), (846) 



where 



d 

c(t) = — \n(a\6c(t)) 



1 d > (85) 

= ---^\a(t)\ 2 + a*a(t). 

These equations show a phase advancement by A, which corresponds to the energy 
separation of states |0) and |1), and both a phase and loss from c(t), which corresponds 
to the phase advance and dephasing from the light. We define 

g c0 (t) = (a\a(t))e tAt pf(t), (86a) 
g 10 (t) = (a\a(t)y At pl°(t), (86b) 

which obey the simpler equations 

ff° = -iga(t)S(t)pl°(t) + (iQ - T)g c0 (t), (87a) 
Q 10 = -iga*S*(t)pf(t). (87b) 

With the phase advances thus removed, the terms g show some damping due to 
spontaneous decay of the atom. When the atom decays, the released photon, phonon, or 
Auger- ionized particle carries away "which-path" information for the qubit, inevitably 
causing decoherence. We have already traced over those lost modes in the derivation of 
the quantum master equation. 

A typical solution of these optical Bloch equations is shown in figure El using 
parameters typical of the 31 P:Si system, as will be discussed in sectional All simulations 
assume F m (t) takes a Gaussian shape with root-mean-square time-pulse-width a P . These 
simulations are performed using a standard adaptive Runge-Kutta solver. It is seen that 
Im{a(£)} rises during the pulse and reaches an asymptotic value in accordance with the 
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Figure 9. (a) Results from a typical simulation of the optical Bloch equations, using 
parameters a = 100, ff\/ K /7 = 20 MHz, 7 = 280 MHz, and r = 300 ns, typical of 
the 31 P:Si system. The real parts of the off-diagonal elements of the density matrix 
are shown, (b) The spectra of the input pulse, cavity response, and 31 P lineshape 
corresponding to this simulation. 



features predicted from both the perturbative analysis of section IP1 and equation (JH2J)- 
The off-diagonal elements g w (t) and g e0 (t) oscillate during the pulse due to the effective 
interaction with the light. 

We may analyze the dynamics of g 10 (t) and g e0 (t) under a similar set of 
approximations under which we approximately solved the first three of the optical Bloch 
equations. We find 

rtVM^vf ^y tf). PS) 

During the pulse, the qubit's off-diagonal element g 10 (t) shows a large oscillation with 
instantaneous frequency of approximately g 2 \aS(t)\ 2 /Q, proportional to the optical 
power inside the cavity. The reduction of the magnitude of this term is due to both the 
divergence of a(t) from a, which is compensated for by returning the factor of (a \ a(t)) 
to Pa°(£), and dephasing due to the loss during the dispersive interaction. It is this effect 
that must be minimized for a high-fidelity gate. 
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Figure 10. A sketch showing how the performance behaves as a function of <& and 
the product of the pulse length and the decay rate, a F T. In the hatched region, the 
pulse length is insufficient to achieve the desired distinguishability. Above that region, 
the fidelity improves somewhat with longer pulses and substantially with higher 



D(t) = e 



|a-aW| 2 /2 



9) 



The total magnitude of the damping to the desired coherence is calculated as 

g 10 (t) 

The final fidelity of our qubit, assuming appropriate single-qubit phase corrections have 
been applied, may be written 

F = l-2\pl°(0)\ 2 (l- lim W). (90) 



t— »oo 

For the calculations presented here, we assume this interaction is being used for 
entanglement distribution, in which case we use p^(0) = 1/2. 

4.. 4-3. Approximate Optimization. We now have enough approximations to estimate 
the optimum values of a, Q, and the pulse width in order to obtain a desired 
distinguishability d = a9 with a maximum fidelity. The roughness of our approximation 
is evident in the loose definition of g 2 \S\ 2 , which we suppose proportional to our estimate 
for the integral over g 2 \S(t)\ 2 , divided by a measure of the length of the pulse, its 
pulse width o~ P . We also assume, quite appropriately, that 9 and L are much less 
than one. For convenience, we use the unitless variable y = Q/T, the pulse detuning 
normalized by the width of the atomic line, which we assume is positive and much larger 
than 1. We also define 

$ 2 r 2 

d l = « rv P $. (9i) 
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We may then summarize our approximations as 

\a\$y 



d 



\ogD ~ — $|a 



l + 2/ 2 + \a\ 2 <S> 2 /Ad 2 m 
1 



(92) 
(93) 



Our goal is to maximize D while holding d constant. 

First, we note that the equation for d may be rewritten as 



where x 



dm _ 
d 

|a|$/2d r 



1 + 



i + (y 



x 



2yx 



(94) 



This equation shows clearly that no matter how large a oc x is 
made, d is always upper-bounded by d m . To reach a desired distinguishability, both $ 
and the pulse width cr P should be sufficiently large to assure d m > d. 

For high fidelity, we must have d m ^> d and work in an unsaturated regime. The 
optimum working condition is to work in the limit where y — > oo and a — > oo, keeping 
constant the ratio 

d 2 



2<- = 

y 

under which condition 
D ~ exp 



\a\$T 



d i 



Ml 



d 2 



2d 2 



(95) 



(96) 



These approximate results are sketched qualitatively in figure [TUJ Quantitatively, 
the critical values of the pulse width and $ to achieve a useful fidelity should be 
calculated numerically, the results of which are shown in the next section. 



5. Calculation of the Dispersive Light-Matter Interaction: Results 

Besides an estimate for the fidelity, the discussion in the previous section provides a 
means of visualizing the output of numerical simulations. The ratio x/y is proportional 
to a/ujQ for systems with a small Purcell effect (either because the mode volume of 
the cavity is high or because u ^> --/). To make this a unitless number, we multiply 
by the normalizing frequency of the simulation. This timescale is the rate at which 
the atom couples to the cavity, g, but we note that in the equations of motion this 
parameter always appears as a product with the ratio yJ~Kpy, which indicates how well 

the cavity is overcoupled. We therefore combine these parameters into g' = gJ k/j; 
this single parameter g' encapsulates all the information about the atom's oscillator 
strength and the cavity geometry. Using this normalization, the parameter ag'/uo, 
indicates the degree to which the interaction is saturated, and therefore we call it the 
saturation parameter. We will plot observed distinguishabilities d = Im a and fidelities 
as a function of this saturation parameter. 

We expect from the analysis above that the distinguishability will rise and reach 
a peak as a function of the saturation parameter. The peak value is d m and occurs at 
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a/u pa 2d m /<&T. Therefore as the pulse length increases, the maximum increases and 
moves to higher values of a/ujQ. If the maximum reached is much larger than the desired 
distinguishability d, then the highest available fidelity is observed; lengthening the pulse 
duration only helps a small amount under this condition. Long pulse lengths do not 
appreciably slow down the proposal for entanglement distribution unless they approach 
the classical communication time between stations, which is about 50 fis for 10 km. 

The behavior of d and F as a function of the saturation behavior should be 
approximately independent of the value of the product aujQ, at least for small values of 
a/u . The product au^ becomes important at high values of the saturation parameter 
due to the Purcell effect. 

In this section, we will see such plots for three important systems. We begin with 
the phosphorous donor impurity in a silicon microcavity, which has an intermediate 
value of $. We will see improved performance in the high-$ system of the flourine 
donor impurity in a ZnSe microcavity. We also consider a typical trapped ion in a 
macroscopic cavity, which may operate in the vicinity of $ ~ 1. 

5.1. 31 P:Si 

The 31 P donor impurity in silicon has become a favored system for several quantum 
information hardware proposals. An early, promising proposal for a quantum computer 
manipulated its electron and nuclear spins using electrostatic gates [2S]. Since then, 
careful measurements of the electron spin decoherence time have been made [2Ej- These 
results indicate that isotopically purified silicon at low temperatures shows extremely 
long electron-spin coherence times approaching 60 ms. 

The spin- 1/2 31 P nucleus is a principal advantage of this impurity for a quantum 
repeater. In a quantum repeater system, very long coherence times are needed, as 
entanglement must be stored for at least the classical communication time over the long 
distances for which such systems are intended (1000-10000 km), and possibly much 
longer depending on the efficiency of the entanglement purification and entanglement 
swapping protocols. A long-distance repeater system will benefit from a many-second 
coherent quantum memory, and probably some form of quantum error correction. In 
semiconductors, only nuclei have shown coherence times this long |2Z|. A critical 
architecture component is therefore the existence of a nuclear spin to which coherence 
may be stored for long periods of time. While proposals for storing electron spin 
coherence in polarized nuclear ensembles have been considered, the degree of nuclear 
polarization required for both suppression of nuclear spin diffusion and high-fidelity 
nuclear memory do not appear to be practical. If a single nucleus such as 31 P in 
Si is used, electron-nuclear transfer may be accomplished by electron-nuclear double- 
resonance techniques, as in the recent demonstration with nitrogen-vacancy centers in 
diamond |28j . Such transfer techniques were first considered in the electron-nuclear 
system of 31 P:Si ^Hj, and have been a strong candidate since some of the earliest 
proposals for experimental quantum information devices [2T?j . 
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Recently, the optical characteristics of this important system have gathered much 
attention. The donor ground state features two Zeeman sublevels, providing the two 
lower levels of the desired A system, and the state |e) is provided by the lowest bound- 
exciton state. It was recently demonstrated that the optical bound-exciton transitions 
are extremely sharp [20], sharp enough to reveal the hyperfme splitting from the 31 P 
nucleus |31j . Such sharp lines allow the measurement [32] of single nuclear spins and 
possibly optical electron-nuclear state transfer, a very rare possibility in semiconductors. 
Unfortunately, this system suffers from silicon's poor optical efficiency The radiative 
lifetime of the phonon-assisted bound-exciton to donor-ground- state transition is about 
2 ms j3S]- The more likely decay channel is Auger recombination, which occurs with 
lifetime 300 ns. This reduces the value of $ for this system by a factor of 10 4 . 
Fortunately, however, microfabrication in silicon is an extremely mature technology. 
The existence of high-quality source material with small optical absorption and the 
wealth of techniques for etching silicon have led to very good lithographically fabricated 
photonic crystal microcavities, with Q-values of 10 6 and mode-volumes on the order of 
(X/n) 3 [HI]. For this system, the Purcell effect makes the system optically active and 
therefore a strong candidate for a CQED quantum repeater. We estimate a coupling 
timescale of g' /2tv = 20 MHz. A Q of 10 6 leads to 7/2vr = 280 MHz. The figure of 
merit $ is then 11. 

Figure [TT] shows the distinguishability d and fidelity F as a function of the saturation 
parameter ol/ujq for this system. For low a/u) , the Purcell effect is not important and 
the dynamics depend little on the product aojQ. The distinguishability and fidelity 
are improved at higher values of a and correspondingly smaller values of u . The 
pulse duration of o~ P = O.lr used in the left plots of figure [TT] is insufficient to allow 
distinguishabilities greater than 1. To increase the maximum distinguishability, the 
pulse must be lengthened at same or higher values of a. The right plots of figure ^TJ 
show how these curves change as the pulse length is increased, each keeping the product 
aw constant at 3162g . Note that the longest pulse considered in this simulation is still 
less than 1 /is. 

For silicon, when the parameters describing the interaction are sufficiently strong 
to allow d ~ 1.6, Auger-limited fidelities greater than 95% are possible. This is sufficient 
fidelity for entanglement distribution, but is likely insufficient for efficient purification 
and entanglement swapping with the measurement-free C-Z gate. For local operations, 
either higher-Q cavities (which have been argued to be theoretically feasible |34j ) or 
another implementation for local quantum logic would assist in efficient purification. 

5.2. 19 F:ZnSe 

The flourine donor in ZnSe is a very promising system. Preliminary measurements [SS] 
indicate that the radiative lifetime of the 440 nm bound-exciton to donor-ground-state 
transition is about 500 ps. Except for the wavelength, this system is therefore optically 
similar to charged quantum dots based on ill-v semiconductors. However, the 19 F:ZnSe 
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Figure 11. The top plots show the distinguishability as a function of the saturation 
parameter for 31 P:Si. In each curve, a and uj are both varied while maintaining a 
constant value of aujo- The bottom plots show the fidelity of the interaction for the 
same parameters. The plots on the left are for a P = O.lr = 30 ns with varying values 
of awo shown. The plots on the right all use acuo = 3162g' but with varying pulse 
lengths listed in terms of the Auger lifetime r = 300 ns. 

system shows some distinct advantages over m-v systems. Bound exciton transitions 
are much more homogeneous in comparison to quantum dots. In comparison to bound 
excitons in GaAs, the higher binding energy and smaller Bohr radius of the ZnSe 
system suggests that it is more robust, easing the isolation of a single impurity in a 
fabricated microstructure. Most importantly, however, the decoherence of electrons in 
GaAs quantum dots or impurities is severely limited by nuclear spin diffusion, as every 
Ga and As nucleus has non-zero nuclear spin. In contrast, the only non-zero nuclear 
spins in ZnSe are the 67 Zn nucleus, which comprises only 4.1% of isotopically natural 
Zn, and 77 Se, which comprises only 7.6% of isotopically natural Se. The 19 F impurity 
is convenient because 19 F is 100% abundant and spin-1/2, like the 31 P nucleus. The 
19 F:ZnSe system is therefore very similar in its nuclear environment to 31 P:Si. 

Perhaps the largest unknown about ZnSe is whether high-Q microcavities may be 
fabricated from this material. Single impurities have been isolated in this material using 
wet-etching techniques [36 j , a promising start toward single impurities in ZnSe photonic 
crystal cavities. Cavities made from distributed Bragg reflectors and microposts are also 
a possibility. We believe it is realistic to expect microcavities with small mode volume 
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Figure 12. On the left are simulation results for F:ZnSe, and on the right are 
simulation results for a trapped ion, both plotted as in figure ITT1 



V ~ (X/n) 3 and moderate Q values of at least 10 3 . 

Using these parameters, simulations of this system yield distinguishabilities and 
fidelities as shown in figure H21 For pulses lasting several nanoseconds and a ~ 100, 
fidelities above 99% are possible. This system is therefore a good candidate for both 
entanglement distribution and the measurement-free deterministic C-Z gate. If the Q 
can be made higher than 10 3 , the fidelity will only improve. 

5.3. Trapped Ions 

Among the longest lived coherences observed in the laboratory are the hyperfine states 
of trapped atoms or ions. Further, a large number of experiments have demonstrated 
the feasibility of quantum logic between multiple ions in a trap, and recent progress 
has been made toward scaling such experiments up to many qubits. For these reasons, 
trapped ions are promising candidates for quantum repeaters. 

Unfortunately, it is very challenging to use a small-mode- volume cavity with an ion 
trap. Therefore, the Purcell factor will inevitably be limited. For typical parameters, 
consider the trapped 40 Ca ion/cavity systems reported in [H3 EH EH] • The atomic levels 
employed in these studies are not exactly the desired A system, but we use values of 
g and 7 from as plausible parameters if such a system were employed. Although 
the cavity Q is quite high (~ 3 x 10 8 ), the large mode volume leads to the estimate 
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$ ~ 1.7. As a result, the fidelity of the system is quite low. Simulation results are 
shown in figure IT21 In this low-Purcell-factor regime the results are nearly independent 
of ac^o, and the distinguishability depends strongly on the pulse- length. However, for 
the desired regime of d ~ 1.6, the fidelity is nearly unchanged after a P > lOr, indicating 
that further elongation of the pulses will allow only a marginal improvement. For the 
ion system, then, internal losses will be as important as external losses in determining 
the final entanglement fidelity. 

6. Discussion & Conclusion 

Most of the components of this proposal for a high-speed repeater are extremely practical 
with current technology. The light source is a filtered, power-stabilized laser with 
pulses that may be generated with normal modulation techniques. The detectors are 
not critical, and do not require abnormally large efficiency or low dark-count rates. 
The required fiber stabilization between repeater stations is available with current 
experimental techniques. 

The only "exotic" element is the effective Hamiltonian of equation (j3J). However, our 
calculations show that this effective Hamiltonian is available in realistic emitter /cavity 
systems as long as the cooperativity parameter $ is large. Such a system is well 
realized by semiconductor microcavity systems in the intermediate coupling regime. The 
19 F:ZnSe system and similar quantum dot systems in photonic crystal microcavities 
with reasonable Q-values show particular promise; here sufficiently large phase shifts 
are available to allow strongly-entangling distinguishabilities with interaction fidelities 
near 99%, resulting in strictly fiber-loss-limited initial entanglement fidelities. Even 
the optically dim 31 P:Si system is reasonable with a sufficiently high-Q microcavity. 
Trapped ion or trapped atom systems may be limited by large cavity mode volumes; in 
such systems the fidelity of the dispersive atom-cavity interaction may only be of order 
80%. Even this system can in principle be well-utilized if sufficiently fast and accurate 
techniques for entanglement purification are employed. 

Although dispersive interactions seem to be strong enough for entanglement 
distribution, a remaining question is whether these systems can be employed for scalable, 
fault-tolerant quantum computation using techniques such as those described in [To] . 
This will ultimately depend on the error- correcting protocols employed. The problem of 
a quantum repeater system is somewhat less stringent since errors of several percent 
during purification and swapping are tolerable ^2]- The use of measurement-free 
C-Z gates with the finite fidelities calculated in the present study should allow some 
rate of purification and swapping to allow scalable long-distance communication. Our 
calculation of communication rates approaching 100 Hz for repeater stations every 
10 km of a total distance of 1280 km was provided as an example; higher speeds 
or longer station-to-station distances are possible. Future theoretical development of 
this proposal will involve optimizing the details of the entanglement purification and 
swapping protocol. 
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Appendix A. Loss Analysis of Deterministic C-Sign Gate 

A wide range of errors are possible for the measurement free C-Z gate described in 
section 13.41 the displacements may be imperfect in magnitude and phase, the single- 
qubit rotations may be imperfect, and of course each interaction of the optical bus with 
a cavity introduces some external loss and some internal loss. For internal loss, the 
fidelity calculations presented in Sections |U and El apply, but in real optical systems 
external losses are always an issue, so we focus on those errors here. 

The first operation in the gate is the dispersive interaction described in section |3] 
To this we must add the displacement of the coherent bus. Although some loss may 
be associated with the displacement, the result of this loss is not very different from 
the result if the loss occurred during the the cavity-light interaction either preceding 
or following the displacement, and therefore we associate loss with the cavity-light 
interaction as described in section 13.21 For example, after the first interaction a loss 
characterized by transmission T\ is assumed to occur; the first displacement in the 
C-Z gate is then altered to \fT[ai{i — 1), where a% would ideally have the same 
phase and magnitude as oiq, but in practice may be slightly different. Since this 
displacement is purely unitary, we may keep track of its operation on the ket resulting 
from equation (fT2j) : 



D[(i - lVTiai 



T^Z^ + ii -IK] . (A.l) 



If the bus state now interacts with a second qubit, the resulting state is 
Tr h U 2 D[(i - 1) Jr^U^piaolUlD^it - l)Jr iai ]U} = 



e 



iTiIm{/3*ai(i-l)} 



T x T 2 p 2 {Z x , Z 2 ) >Q 2 (Qi(p))< ^JT 1 T 2 /3 2 (Z 1 , Z 2 



(A.2) 



where we have traced over lost photons. The superoperator Q 2 (p) is substantially more 
complicated than Q\{p). In general, at the nth step in the gate, we may write Q n (p) m 
the Z-basis as 

Qn{p) = \jk)(jk\p\mn)(mn\ x 



jk,mn 

exp 



, , (A.3) 



1 - T n )T n . lTn . 2 ■■■[ imk) f n{m > n)l + ilm{m k)P n (m, n)} 



where j and k enumerate the eigenstates of Z\ and Z 2 ^ as do m and n. It is worthwhile 
to point out here that Q n will in general contain state-dependent phase shifts, and that 
these will contribute to the C-Z gate as well as the Berry phases which accrue during 
displacements. 
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This type of analysis may be carried through to the end of the gate. The full 
sequence of operations is 

U 4 D[(1 - i)^TiT 2 T 3 a 3 ]U 3 D[(-l - i)^I\T 2 a 2 ]U 2 D[{-l + z)/^]^, (A.4) 

where U3 and U4 correspond to couplings back to qubits 1 and 2, respectively, but may 
have different transmissions Tj and angles 9j due to imperfect coupling in the optical 
circuit. As a result of this sequence, the total state-dependent Berry phase that develops 
is given by 

lm{(i - 1)^^(^)01 - (i + \)T{T 2 p* 2 {Z x , Z 2 )a 2 - (i - 1)T X T 2 T^{Z X , Z 2 )a 3 }. 

At the end of the gate, the coherent bus photons are lost; i.e we may use T 4 = with 
equation (|A.3j) to find the final superoperator Q±{p) after losing the bus state. 

With this formalism, several errors in the C-Z gate may be analyzed in general. 
Here we suppose that the displacements are accomplished perfectly (every aij is the same 
in magnitude and phase), and every interaction with every cavity is exactly the same, 
hence studying only errors due to optical loss. We call this the semi-ideal case. Also, 
we work to second order in 9, as this is the lowest non-vanishing order after correction 
of single-qubit phase shifts. 

In this semi-ideal case the coherent state amplitudes witnessed by the bus state 
during the gate are 



(3x{Z x ,Z 2 ) = In | 2 1 

HZi,z 2 ) 

^(Zt,Z 2 ) 

fa{Zi,Z 2 ) 

and the Berry phase simplifies to 

9 





2 


a 






2 


a 






2 


\a 






2 


a 





%—Z\ 

2 8. 



[l+i + 2Z 1 Z 2 ) 



-Z 2 -j[l + (l+i)(l + 2Z 1 Z 2 )} 
j(l+i){l + Z 1 Z 2 



(A.60) 
(A.66) 
(A.6c) 
(A.6d) 



\V a \ 



|[1 +T + T 2 ] + -[(i+T)(Z 1 +TZ 2 )] + -[3T 2 - 1 + 2T(1 + 2T)Z 1 Z 2 ]\.(A.7) 

We see that this may be immediately decomposed into a global phase (due to constant 
terms), single-qubit rotations (single Z terms), and a nonlinear term going as Z\Z 2 . 
This last term is part of the nonlinearity leading to the C-Z gate. 

However, extra phases occur due to the superoperators Q n , which resulted from 
"which-path" information carried away by lost photons or the optical bus. To second 
order in 9 in the semi-ideal case, 



(A.8) 



Q*(Q 3 (Q2(Qi{p)))) = 

e i|a| 2 [(l-T)(Zi+TZ 2 )0/2+T 2 (l-2T)ZiZ20V4]^^ e -i|a| 2 [(l-r)(Zi+rZ 2 )e/2+T 2 (l-2T)ZiZ 2 eV^ 

where we have seperated out pure rotations and C-Z terms from a distortion operator 
D(p), which may in general be generated from equation (|A.3)1 in the computational 
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basis as D(p) = Y,jk,mn'Djk,rnn\jk)(jk\ p\mn)(mn\. The matrix V is Hermitian with 
unity diagonal. In the semi-ideal case, to lowest order in 9, the off-diagonal elements 
are 

V ++ , + _ = P__,_ + = e -[ni-TV^i-T)]MlV2-M 2 | 2 /4 ; (A-9a) 
= V__^_ = e-I (i-T 3 )-Hr(i-r)]|atf| 3 /2-|atf»IV4 > (A . 96) 
!?++,__ = !?+_,_+ = e -( 1 - T2 )( 1+T )l^l 2 /2. ( A. 9c ) 

Terms of order 9 A are present to show that there is distortion even if there is no loss in 
the system (T = 1). 

For the operation of the C-Z gate, we presume the single-qubit rotations from the 
semi-ideal case expi(l + T 2 )(Z\ + rfZ 2 )\a\ 2 9 /2 are removed. Imperfect correction of 
these phase shifts may be an important form of error, as these shifts are of order \a\ 2 9, 
which will be much larger than tt. Neglecting these shifts and the global phase shift, 
the resulting unitary operation from this gate is the very simple 

U C -z = ex P (iT(l + T)^Z 1 Z 2 y (A.10) 

To make this a C-Z gate, we must choose T(l + T)|a#| 2 = it, and add single-qubit 
rotations exp[— i(n/4:)(Z 1 + Z 2 )}. 

For a general case, the distortions due to errors are more computationally convenient 
to express as a sum of Kraus operators. Let A m be the eigenvalues of the matrix X>/4, 
and define A = H 2 T, where H 2 is the usual 2-qubit Hadamard transformation and the 
columns of T are the orthonormal eigenvectors of T>. Then the distortion superoperator 
may be written 

D(p) = J2D mP Dl (A.ll) 

m 

where 

D rn = \J X m (A ++>m + v4_| — )m Z 2 + A — \. >m Z\ ). (A.12) 

For example, we used such a description for the Q\ superoperator, with only two noise 
components, D + = ex P ^i£i/2 and L>_ = yf\ZZ\ expiZi£i/2. The eigenvalues 
X m are in general calculated numerically. The only operator with an identity element is 
Do, corresponding to the largest eigenvalue Ao- Therefore the deviation from unity of 
Ao characterizes the fidelity of the gate. For 9 2 <^ vT, which is the appropriate regime 
for the systems considered in this paper, the largest eigenvalue A is approximately 
(1 + exp(— 7ta/T/4)) 2 /4. Numerical results for several values of 9 are shown in figure[3 
The distortion operators D m commute with the desired C-Z gate. To generate a 
controlled-NOT, or C-X gate, we must simply rotate the target qubit, qubit 2, with 
exp(— iirY/4). Therefore the final noisy operation describing a C-X gate in the presence 
of optical loss is 

^c-xEVM-x, (A.13) 

m 
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+ A.^ m Z l X 2 ). (A.14) 

These operators fully describe the error model which may be used in an analysis of 
entanglement purification schemes or error- correct ion techniques. 

Appendix B. Phase-Shift Calculation with Finite-Order Cumulant 
Expansion 

In this Appendix we detail an approach for calculating the state-dependent phase shift 
due to a single atom-cavity interaction in the dispersive limit. Here we calculate the 
phase shift and internal loss assuming the atom begins in state |1). 
We perform a finite-order cumulant expansion, in which we write 

x(v, t) = x(v, 0) exp[^(e, 77, t)}. (B.l) 

The function \I/(e, r], t) vanishes at e = and t = 0; for e > it may be expanded in a 
Taylor series 

e 2 

*(e, 77, t) = ettifa, t) + -y 2 ( V , t) + ... (B.2) 
This may be compared to the Taylor series for x(v^)' 



X{v,t) = Xo{v,t) 



e 2 

l + e^! + -(^ 2 + ^ 2 ) + ... 



(B.3) 



the latter of which we actually calculate to finite order. 

To find each order, we solve the Liouville-von Neumann equation for the atom- 
cavity coupling. The iterative solution is 

00 

p(t) = £ (B.4) 

k=0 

where po — p(t — 0)- We again work with Laplace transforms, with which we find 

Pk{s) = — f 7T- f W~ + -£[Pk-i(s)], B.5 

s J 2m J 2m s — q\ — q 2 s 

where in the inverse transform we must close the s contour to the right of all qi,q2- 

Using equation (|37cjl . pk(s) may be found from 



Pk{s) = ig^J2f 



dq 
2mL 



Fm(q) + , KM t - Pk-i(s-q) 

1 a ™° 1 Z a m a > 

q + z q + z* s 



+ - s C\p k - X {s)\, (B.6) 



where for compact notation we define z = iujq + 7/2 and F Q (q) = — The first few 
terms are 

Po(s) = - \oc)(a\ ® o~o + 



s 



Pi(s) = iqJTii a—-^- r \a)(a\ a + — h.c. 

r w I s(s + z) 1 A 1 

p 2 ( s ) = - g\^AW(s)[ala-, a \a)(a\ a + ] + h.c.} - ^ 



t Pi{s) 
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-^rH[^m(s)[a^a-,a|a)(a|(T + ] +h.c.| - 



2 

2— Re{4 N (s)}|a| 2 \a)(a\a 2 

ST 



F m (q)[A m (s -q) + A* N (s - q))a* \a)(a\ a m 
+ F m (q)[A m (s - q)a\ n a \a)(a\ + A* m (s - q) \a)(a\ a*a m ] \aa' 



where 

A®(8) = 1 ^- F MF m {s-q) 

ml J / 2iris{s-qy{q + z*){s-q + zy { ' 

The two lowest non- vanishing terms of the Laplace transform of \I/(e, rj, t) are then 
found as 

i* 2 fa, s) = Xo" 1 ^) I^AjurfaWs)} = +<7 2 «A ( N } (^*a - c.c. (B.8a) 

~¥ 8 fo, s) = Xo "(^) Tr{^ OUT (77)p 3 ( S )} = -^^W^a - c.c., (B.86) 

where "c.c." refers to complex conjugate. Again, we are interested only in the asymptotic 
t — > oo limit. All poles of A(s) with negative real part may be ignored. Poles with 
imaginary part also damp out since the pulse has finite duration. Only the s = pole is 
important. We write q = iu and note that the denominator (s — iuj)~ l may be written 
i/uj + ttS(u), resulting in 



2tt (jj\z — iuor 2 \zr 



From these expressions we obtain the estimates discussed in section 14.31 
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